Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



SIMULATION OF SOFT X-RAY EMISSION LINES FROM THE MISSING BARYONS 

Taotao Fang^ *^, Rupert A. C. Croft\ Wilton T. Sanders^, John Houck^, Romeel Dave'', Neal Katz"', 

David H. Weinberg", and Lars Hernquist' 

(Received; Revised; Accepted) 
Draft version February 2, 2008 

ABSTRACT 

We study the soft X-ray emission (0.1 - 1 keV) from the Warm-Hot Intergalactic Medium (WHIM) 
in a hydrodynamic simulation of a Cold Dark Matter universe. Our main goal is to investigate how 
such emission can be explored with a combination of imaging and spectroscopy, and to motivate future 
X-ray missions. We first present high resolution images of the X-ray emission in several energy bands, 
in which emission from different ion species dominates. We pick three different areas to study the 
high resolution spectra of X-rays from the warm-hot IGM: (a) a galaxy group; (b) a filament and (c) 
an underluminous region. By taking into account the background X-ray emission from AGNs and 
foreground emission from the Galaxy, we compute composite X-ray spectra of the selected regions. 
We briefiy investigate angular clustering of the soft-X-ray emission, finding a strong signal. Most 
interestingly, the combination of high spectral resolution and angular information allows us to map 
the emission from the WHIM in 3 dimensions. We cross-correlate the positions of galaxies in the 
simulation with this redshift map of emission and detect the presence of six different ion species (Ne 
IX, Fe XVII, O VII, O VIII, N VII, C VI) in the large-scale structure traced by the galaxies. Finally 
we show how such emission can be detected and studied with future X-ray satellites, with particular 
attention to a proposed mission, the Missing Baryon Explorer, or MBE. We present simulated 
observations of the WHIM gas with MBE. 

Subject headings: cosmology: theory — diffuse radiation — large-scale structure of universe — X-rays: 
diffuse background 



1. INTRODUCTION 

Cosmological hydrodynamic simulations provide us 
with quantitative predictions for the state of baryonic 
matter in the u niverse (see, e.g, ICen fc Ostrike3ll999l: 
iDave et al.l200H) . At high redshift (^ ^ 3), most baryons 
are diffuse, with a density close to the cosmic mean, and 
are located in the intergalactic medium (IGM). They 
can be probed readily using Lya absorption systems 
in the optical sp e ctra of high redshift quasars (see 
e.g., iCen et al.l 1X9941: iZhan g. Ann inos. fc NormanI 
19951 i Hernauist. K atz. Weinberg, fc Jordil Il996t 
Weinberg e t al. 1997). At low redshift, the process 



of gravitational collapse has caused a minority of 
baryons to condense into structures such as stars, 
galaxies, groups and clusters of galaxies. While a sub- 
stantial fraction (~ 30%) of baryons are trapped within 
the gravitational potential wells of these large scale 
structures in the form of hot intracluster or intragroup 
gas, the remaining baryons are left in intergalactic space 
and are predicted to form filamentary structures seen in 
cosmological hydrodynamic simulations. Some of these 
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structures contain baryons with relatively low tempera- 
tures (T < lO'' K), and therefore enough residual neutral 
hydrogen to be d etectable as the low-redshift Lya forest 
(see, e.g.. iDav^. Hernauist. Katz. & Weinberg 199^ 
iPenton. Shull. fc St ockc 2000). However, a larger frac- 
tion have been shock-heated to temperatures between 
10^ — 10^ K. Given these temperatures, these baryons, 
which form the so-called "warm-hot intergalactic 
medium," or "WHIM," are best probed using UV/X-ray 
observations. 

Recently, important progress has been made in 
the detection of absorption features produced by 
highly ionized metals in this WHIM gas in the 
UV/X-ray spectra of background quasars. In the 
UV band, a significant number of O VI absorption 
lines have been seen with the Far Ultraviolet Spec- 
troscopic Explorer (FUSE ) and the Hubble Space 
Telescope (HST) (see, e .g., 'Savage. TrioD. fc Lu ^199^ 
|T ripp. Sav a ge, fc Jenkin ^ F oOO: TriDD & Savage 20o3 
ISimcoe. Sargent, fc Rauch "1200^1 . The distribution 
and derived properties of these O VI lines are con- 
sistent with p redictio ns from simulations (see, e.g.. 
Fang fc Bryan* '2001'; 'Cen. TrioD. ©striker, fc JenkinJ 
' 01; Chen. Weinberg. Katz. fc Dave 2003). In the X- 
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ray ba nd, IFang et al. ( 20021. .Mathur. Weinberg, fc Cheil 
( 2002[l . iCagnoni ' (,2002^) and IMcKernan et alJ 
( 2003}) reported on the detection of interven- 
ing O VII and/or O VIII absorpti on lines with 
Chand ra and XMM-Newton. iNicastro et alJ 

f2002^, 'Fang, Sembach, & Canizar^ ||2003D and 



R,asmussen. Kahn. fc Paerels (200,1) also reported 
the detection of z w X-ray absorption lines. These low 
redshift lines may be attributable, at least in part, to 
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the WHIM gas in our Local Group. 

The absorption hne method does have the advantage 
of directly probing the absorbing ions, as the detected 
line strength is proportional to the column density and 
so to the ion number density. However, it suffers the 
limitation of being only able to probe one dimensional 
information, along limited lines of sight, and detectabil- 
ity is also constrained by the flux of background sources. 
In order to fully reveal the three dimensional structure 
and physical properties of the WHIM gas, we need to 
study its emission with imaging/spectroscopic methods. 
In this paper, we will focus on the X-ray emission from 
the WHIM gas. In general, X-ray emission from a hot, 
diffuse plasma contains two parts, a continuum which 
is contributed by various emission mechanisms, and line 
emission from metals in the IGM. While we will briefly 
discuss the continuum emission from the WHIM, this pa- 
per will be devoted mostly to studying the line emission. 

Several authors have extensively studied the broad 
band soft X-ray emission (0.5 ~ 2 keV) from the 
diffuse background, particularly the X-ray emission 
from the WHIM gas, using hydrodynamical simu- 
lations (see, e.g.. iCen Kang. Ostriker. fc RviJ 119951 : 
IPhillips. Ostriker. fc Cenll2nniHCroft et al.ll2n0ir) . Much 
effort has been put into investigating the overall intensity 
of the emission, and how to distinguish it from the two 
dominant sources of the soft X-ray background: X-ray 
emission from active galactic nu clei (AGNs) and from the 
hot gas in our Galaxy (see. e.g.llCiintz fc Snowdenll2fl0nl: 
iKuntz. Snowden. fc Mushotzkvll200lD . Simulations pre- 
dict that the mean intensity of emission from the WHIM 
gas is between (2 - 4) xlO~^^ ergs s~^cm~^deg~^, only 
about 5 - 15% of the total extragalactic X-ray emis- 
sion between 0.5 - 2 keV l|PhilliDS. Ostriker. fc Cenl2001l: 
iCroft et all 12001(1 . While it seems to be a difficult 
task to detect most of this WHIM emission with cur- 
rent X-ray telescopes, it appears feasible to detect the 
high intensity tail of the WHIM X-r a y inte nsity dis- 
tribution. For instance, IScharf et al.l l|200ClD reported 
on the presence of a possible filament with intensity 
» 6 X 10"^^ ergs s~^cm-^deg ~'^ with the ROSAT PSPC; 
and iMarkevitch et all l|2003li and iKaastra et all lj2003l) 
reported detecting O VII and/or O VIII from an ex- 
tended, diffuse structure , possibly from WHIM gas; and 
iMcCammon et aTl ||2002D recently performed a large-field 
(~ 1 sr), high spectral resolution (5 - 12 eV) observation 
with a sounding rocket and detected C V/0 VII/0 VIII 
lines from the diffuse background, some fraction of which 
may be from the WHIM at zero redshift. 

We expect that future X-ray missions will enable 
high resolution imaging/spectroscopic observations of 
the emission from the WHIM gas. We explore such a 
scenario with a high resolution hydrodynamic simula- 
tion in this paper. Our purpose is to address two im- 
portant questions: with high spectral/spatial resolution, 
what does the X-ray emission from the diffuse, hot gas in 
intergalactic space look like; and are proposed future X- 
ray missions capable of conducting useful observations? 
Particularly, we will focus on a proposed mission called 
the Missing Baryon Explorer, or MBE. Given its large 
field-of-view (29.5' x 29.5'), high spectral (- 4 eV) and 
moderate angular 5' x 5') resolution, we find that 
MBE will be capable of revealing the X-ray emission 
from the WHIM gas. In addition, we will explore the 



clustering properties of the WHIM gas, and its correla- 
tion with ne arby large scale st r ucture s traced by galaxies. 

Recently, lYoshikawa et al.l l|2003f) examined the de- 
tectability of the X-ray emission from the WHIM gas 
independently with a similar numerical approach. They 
presented a detailed study of O VII and O VIII emis- 
sion lines from WHIM. They also investigated a dedi- 
cated X-ray mission, DIOS (Diffuse Intergalactic Oxy- 
gen Surveyor). It turns out that DIOS and MBE very 
similar properties. While there are some differences be- 
tween our work and th eirs, such as the selection of the 
energy bands ifYoshikaw a et al.ll200,'^ concentrated on 0.5 

- 0.7 keV while we look at a broader band between 0.1 

- 1 keV) and the choice of the IGM metallicity models, 
where they can be compared our results agree with each 
other reasonably well. 

Our plan for this paper is as follows. We give a brief 
explanation of the simulation method and our treatment 
of the simulation data in A full study of the emis- 
sion from the WHIM gas is presented in ^ followed by 
measurement of its auto- and cross-correlation proper- 
ties, making use of the galaxy distribution in the simu- 
lation for the latter. We investigate the detectability of 
emission with the proposed MBE X-ray telescope in fJS] 
and [jSjis our discussion and summary. 

2. SIMULATION 

We use a cosmolqgical hydrodynamic simulation stud- 
ied by iDave et all l|20nil ^ (their model Dl). It was run 
with PTreeSPH, or parallel tree smoothed particle hydro- 
dynamics (see Dave, Dubinski. & Hcrnquist 1997 ). We 
refer t he reader to these two papers as well as Gro ft et alJ 
(|2001|) for more details. The model simulated is a Gold 
Dark Matter (CDM) universe with a cosmological con- 
stant (ACDM), with cosmological parameters JIa — 0.6 
and il,n = 0.4. The parameter h = 0.65, where the Hub- 
ble constant is Hq — 100ft. km s~^Mpc~^. The simulation 
box size is 50/i^^Mpc with a 7ft.^^kpc spatial resolution 
(comoving units, equivalent Plummer softening). The 
baryonic mass resolution is 8.5 x 10^ M©, where Mq is 
the solar mass. Cooling from H and He was included, as 
w ell as star fo r matio n and mild feedback. As discussed 
in IDave et al.l l|2001[ ). the feedback consists of thermal 
energy which is released into star forming regions, which 
are of high density. The energy is rapidly radiated away 
and has little effect on the IGM. The simulation was run 
from redshift z = 49 to the present. We make use of 
27 of the outputs, roughly logarithmically spaced in red- 
shift, ranging from z = 5.98 to z — 0. 

For each gas particle, we calculate the X-ray emissiv- 
ity ba sed on a Raymond-Smith code ( Raymond & Smit^ 
Il977r . Note that we have first recalculated the SPH 
densities for each X-ray emitting particle after exclud- 
ing the cold par ticles from the density estimator (see 
ICroft et 311120011 for details). This decou pling of hot and 
cold phases (see also iPearce et al.l|200nl for example) is 
necessary in order avoid overestimating the densities of 
hot particles which are clo se to cold particles (see also 
iSpringel fc Hernauisdl200^ . 

We adopt a density-dependent metallicity, Z (x ^ p/'p, 
to mimic the d e nsity- m etallicity relation obta i ned b y 
IGen fc Ostrike? (Tg^, ISnringel fc HernauisH ll2f)?)l . 
and Aguirre et al. (2001a,b,c) in their simulations. Here 
p and p are the density and the mean density of the uni- 
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verse. We normalize this relation by setting Z = 0.005^0 
at p = p to match Lya forest measurements, and limiting 
the metallicity to be Z < O.BZq in much more overdense 
regions to match the observational results from galaxy 
clusters in the local universe. Here Zq is the solar abun- 
dance. The Raymond-Smith code includes eleven heavy 
elements, namely C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, 
and Ni. Table 1 lists the wavelengths and the corre- 
sponding energies of the strongest emission lines we are 
probing with our simulation. Assuming coUisional ion- 
ization equilibriu m, Figure [T] shows the em issivity for 
each ion species l)R,avmond fc Smithlll977|) . where i rep- 
resents ion species Xi. We give the strongest line (the 
resonance line in the He-like triplet), and for Fe XVII we 
select the 825.79 eV hne. 

In Figure 12 we display three template X-ray emission 
spectra for hot gas at different temperatures: (a) T = 10^ 
K; (b) T = 10^^-^ K; (c) T = 10^ K. We use a metallic- 
ity of 0.1 Zq, which is typical for gas w ith overdensities 
5 a few hundred based on simulations (* Cen & Ostrikerl 
((19991 Sorineel fc He rnauist ( 2003)). Here b is defined 
as (5 = jo/p — 1. Figure HJi and 12) show that at low 
temperatures the spectra are dominated by C, N and 
O VII emission lines. At high temperatures O VIII and 
Fe emission start to take over. At temperatures above 
10^ K, the spectra are completely dominated by Fe and 
Ne emission lines. 

Table 1: Emission Line List 
Ion Wavelength (A) Energy (eV) Line^ Reference^ 



C V 


41.47 


298.96 


f 


a 




40.73 


304.41 


i 


a 




40.27 


307.89 


r 


a 


C VI 


33.73 


367.54 




b 


N VI 


29.53 


419.79 


f 


a 




29.08 


426.30 


i 


a 




28.78 


430.80 


r 


a 


N VII 


24.78 


500.35 




b 


VII 


22.10 


560.98 


f 


a 




21.80 


568.62 


i 


a 




21.60 


573.95 


r 


a 


VIII 


18.97 


654.00 




b 


Fe XVII 


17.10 


725.22 









17.05 


727.14 




c 




16.78 


738.88 




c 




15.01 


825.79 




c 


Ne IX 


13.70 


905.10 


f 


a 




13.55 


915.03 


i 


a 




13.45 


922.00 


r 


a 
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Fig. 1. — X-ray emissivity for eight ion species. Among them, 
O VII, O VIII, and Fe XVII show the highest peak emissiv- 
ities between 10^ — 10^ K. Emission models were taken from 
■Raymond fc SmithI 



combination dominates at E > 0.4 keV, while thermal 
bremsstrahlung radiation dominates over the entire X- 
ray band when T ~ 10^ K. The total continuum emissiv- 
ity e^, as a function of photon e nergy E', can be written 
as l|Mewe k. Gronenschildlll98H) : 



t,{E) oc G,{E,T)ET 



-1/2 2 



exp{~E/kT), (1) 



Three types of continuum emission make contributions 
to the total X-ray continuum we discuss here: thermal 
bremsstrahlung radiation, or free-free emission; radiative 
recombination radiation, or free-bound emission; an d 
two-photon decay (see. e.g., lR,avmond fc Smitihlll977t) . 
Their relative contributions vary with temperature and 
photon energy. For instance, at T ~ 10^ K radiative re- 



where T is the gas temperature, Ue is the electron den- 
sity, k is Boltzmann's constant, and Gc{E,T) is the to- 
tal effective Gaunt factor, which includes all three con- 
tinuum emission processes (jGronenschild fc Mewdl 978). 
By integrating ec{E) over energy E, the integrated X- 
ray continuum emissivity is oc niT", where a = 1/2 for 
high temperature gas, in which bremsstrahlung radiation 
dominates, and a = —3/2 for gas with T < a few xlO^ 
K. 

3. SOFT X-RAY SPECTRA OF THE WARM-HOT 
INTERGALACTIC MEDIUM 

3.1. The IGM: simulated maps and spectra 

In order to go from the three dimensional particle dis- 
tributions in the simulation outputs to sky maps and 
spectra, we rely on the technique o f stackin g simulation 
boxes along the hghtcone (see e.g., Ida Silva et al|l2000t 
ISpringel. White, fc Hernamsil l200il IGroft et all I2001D . 
We place simulation boxes one behind the other until 
we reach the comoving distance corresponding to a given 
redshift (z = 1 in the first case). The output file cor- 
responding to the closest redshift is used at each point 
(we use 11 different outputs between z = and z = 1). 
The boxes are ran domlv rotated, reflected and translated 
;rowijj(g^faW jti ^200H) so that periodic structure does not re- 
peat along the line of sight. Because of this, only struc- 
ture on scales smaller than the box is preserved, and we 
only consider structure on these scales when measuring 
clustering statistics. For example, the median redshift of 
emission is z = 0.45 in the 0.5 — 2 keV band, so that we 
cannot accurately simulate structure with angular scales 
> 100 arcmins. 

The X-ray emissivity in a range of 2eV wide bins in 
energy is computed for each particle, and these quanti- 
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Fig. 2.— X-ray emissivity for hot gas at three temperatures typical of the z = IGM: (a ) T = lO'' K; (b) T = lO'^'^ K; (c) T = lO'^ K. 
The gas has a metaUicity of 0.1 Zq. Emission models were taken from lRavmond fc Smitlj 119771) . 



ties are allocated to a grid using the projection of the 
SPH kernel on the plane of the sky (see below for more 
details). The grid is a three dimensional datacube, with 
two angular axes and one axis for photon energy. 

We first examine the angular distribution of the emis- 
sion in a relatively broad spectral band. In Figure 01 we 
show a map of a 200' x 200' simulated region. The area 
has been divided into a 512x512 grid, so that each pixel 
has side length ~ 23". The X-ray emission includes both 
continuum and metal line emission from the hot IGM, in- 
tegrated between 0.1 - 1 keV.ICrpft ct al. ( 2001) showed 
that the majority of the X-ray emission from the IGM re- 
ceived at z = co mes from gas with z < 1.0 (about 90%, 
see Figure 11a of lGroft et al.ll200'l|) . Due to dilution by 
projection, filaments that would be apparent in three di- 
mensions do not appear readily here. The intensity varies 
from 10~^ — 10~^ erg s~^cm~^deg~^, in patches of sky 
dominated by emission from the hot IGM in groups or 
clusters of galaxies, down to 10^^^ erg s~"'^cm~^deg~^, 
in regions of low density. We label three (5' x 5') ar- 
eas that we will study in detail spectroscopically in the 
following sections: area (a), with an average surface 
brightness of lO"'' erg s~ cm~^deg~^, representing X- 
ray emission from a galaxy group; area (b), with a mean 
surface brightness of 10~^ erg s~^cm~^deg~^, represent- 
ing X-ray emission from a superposition of filamentary 
structures; and area (c), with a surface brightness of 
10"^"^ erg s~^cm~^deg~^, an underdense and void-like 
region. 

Figure 01 shows the soft X-ray emission from the same 
sky area but in four narrower energy bands. Since the 
X-ray continuum emissivity exponentially decreases with 
energy, in most energy intervals the flux from the X-ray 



emission lines dominates over the continuum emission. In 
panel (a) emission lines from carbon, particularly C VI 
and C V triplet dominate from 250 to 350 eV. Here we 
assume that most X-ray emission l ines come from the 
low-redshift IGM. iGroft et all l|2001() showed that X-ray 
emission from the IGM peaks at z ~ 0.2 and decreases 
slowly to high redshift. Panel (b) shows the band (350 
- 425 eV) in which nitrogen lines dominate. Panel (c) 
is for the band (550 ~ 600 eV) for which oxygen lines 
dominate (mainly O VII triplet). Panel (d), the highest 
energy is in the band where emission from iron is very 
strong, mainly the Fe L complex. It can be clearly seen 
from the emission maps in these four narrower energy 
bands that low energy emission is rather more diffuse 
(panels a and b) than high energy emission, as would be 
expected from a superposition of filamentary structures. 
On the other hand, at high energy emission comes from 
more obviously discrete sources, such as collapsed regions 
which form groups of galaxies. 

To study the spectroscopic properties of the soft X-ray 
emission from the WHIM gas, we divide the simulated 
patch of sky into 40 x 40 grid cells, with each cell ex- 
tending over an area of 5' x 5' ^. In every cell, we obtain 
the X-ray spectrum of each gas particle by applying the 
Raymond-Smith code, based on its temperature, den- 
sity, metallicity and redshift. The peculiar velocities of 
the gas particles are used in order to produce a redshift 
space spectrum. We then combine all the spectra to ob- 
tain a synthesized spectrum in each cell, as mentioned 
above. 

^ The pixel area of 5' X 5' is the spatial resolution of the Missing 
Baryon Explorer (MBE), a proposed X-ray mission which we will 
discuss later. 
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Fig. 3. — The simulated 200' X 200' region. The whole region is divided into 512x512 grid cells, and the X-ray intensity is each cell is 
obtained by integrating the X-ray flux between 0.1 — 1 keV in that cell. We select three regions (shown as white boxes) to study three 
typical types of X-ray emission: (a) hot gas in a galaxy group; (b) diffuse X-ray emission in a filament; and (c) a void. 



Figure El shows 10 random samples of such spectra 
picked from the brightest 10 % of pixels over the range 
0.3-0.9 keV. A large number of emission lines are visible 
in many of the panels. Having information from the un- 
derlying simulation allows us to identify lines easily. Our 
scheme for doing this automatically is as follows. When 
creating the spectra, for each energy bin, we keep track of 
the redshift of the particle that was the largest contribu- 
tor of X-ray flux to that bin. When plotting the spectra, 
we divide the energy of bin z by 1 -|- , where Zi is the red- 
shift of the maximum flux particle. If the resulting value 
is within 5 eV of one of the known ion rest wavelengths 
(see Table 1), and the pixel corresponds to a peak in the 
spectrum, we label it with the name of ion species, and 
the redshift Zi. It can be seen that many of the features 
come from the same redshift, but are signposted by differ- 
ent species. For example, the top left panel shows that 



a number of lines (N VII, O VIII, Fe XVII and Ne IX) 
arise from hot gas at redshift ^ 0.35. These types of 
spectra can potentially provide us with important infor- 
mation that could help us identify emission lines with ion 
species and redshifts in real observations. 

3.2. Contribution from the Galactic Foreground and 

AGNs 

ROSAT llSnowden et alJ Il995l 119971) and other 
X-ray missions (se e, e.g., iMcCammon et al.l l2002t 
iSanders et alJ 120011) show that the soft X-ray back- 
ground in the band 0.1 - 1 keV is comprised of at 
least three components: (1) unabsorbed X-ray emission, 
referred to as the "local hot bubble" (LHB), f rom a 
cavity in t he local interstellar medium (.SnowdenlllQQSt 
iSfeir et al.lil999,) : (2) hotter absorbed X - rays, referred to 
as "transabsorption emission" (TAE, K untz fc Snowdei] 
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FlG. 4.— The simulated 200' X 200' region in four energy bands: (a) 250 - 350 eV; (b) 350 - 425 eV; (c) 550 - 600 eV; (d) 700 - 800 eV. 
Emission lines from C, N, O, and Fe dominate in each band, respectively. 



l2000HSnowden. Frevberg. Kuntz. fc Sand ers"2'000^ . from 
gas in either the Galactic h alo or the intragroup 
medi u m of the Local Group l|Rasmussen fc PedersenI 
|2P07, 'Rasmussen. Kahn. fc Paerels"200 3D, and (3) an 
extragalactic X-ray background ( Has inger et al.l Il998t 
iMushotzkvTGowie. Barger. fc ArnaudI 1200(1) Simulta- 
neous fits to ASCA and ROSAT data arrive at the same 
components l|Mivaii et alJll998|) . 

It is well known that the extragalactic X-ray back- 
ground is dominated by discrete sources such as active 
galactic nuclei (AGNs). ROSAT deep surveys resolved 
~70 - 80% of the soft X-ray background (0.5 ~ 2 keV) 
into discrete sources at a flux limit of 10~^^ erg cm^^s^^ 
ijHasinger et alJ Il998() . The Chandra deep survey, 
which had a flux limit of 3 x 10~^^ erg cm~^s~^, 
resolved an additional 6 - 13% int o point sources 
l|Mushotzkv. Cowie. Barger. fc ArnaudI l200a) . This 



leaves room for ^--^5 - 25% to come from diffuse emission 
from hot gas in the IGM, which is what we concentrate 
on in this paper. 

To obtain the composite spectrum of the Galactic fore- 
ground and the AGN background, we use the following 
simple models of the components. (1) For the LHB, we 
adopt a thermal plasma model with kT = 0.1 keV and 
emission measure EM = 0.019 cm^^ pc. The metal 
abundance is taken to be solar, except for Mg, Si, and 
Fe, which have an abundance of 0.28 Zq. (2) For the 
TAE, we adopt a thermal plasma model with kT = 0.227 
keV and emission measure EM = 0.003 cm^^ pc. This 
model has solar abundances and is absorbed by a hy- 
drogen column density of Nh = 1.8 x 10^° cm~^. The 
parameters for these two thermal components were ob- 
tained fro m fits to the calorimeter s ounding rocket data 
reported in lMcGammon et all l|2002f) using the MEKAL 
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Fig. 5. — 10 random sample spectra taken from the brightest 10% of pixels. Note that the y-axis scale changes between each plot. We 
have plotted a gray bar at count level 10~* photons cm~^ pixel"^ keV~^ so that the panels can more easily be compared. The pixel 
size is 5' X 5' . The redshifts of material contributing to emission lines are shown 



spectral models in XSPEC (3) The third component 
represents the AGN. We used a spectrum constructed by 
D. McCammon (private communication) to reflect all the 
known ROSAT and Chandra data on the resolved AGN 
con tribution of the soft X- ray background (see Figure 
15, iMcCammon et aTl 120021) . In Figure 6 we show the 
composite spectrum of these three components. 

3.3. X-ray spectra of Characteristic Regions 



see |http://heas arc. gsfc.nasa.gov/docs/xanadu/xspec/ 



We select three regions (shown as white boxes in Fig- 
ure|21l to study three typical types of X-ray emission: (a) 
hot gas in a galaxy group (Figure E^); (b) diffuse X-ray 
emission in a filament (Figure Eb); (c) a void (Fig- 
ure[3;). In each panel we show three curves: the red line 
is the emission from the IGM, green is the composite 
spectrum of AGNs and the Galactic foreground, and the 
black line represents the total emission. We also identify 
the most significant IGM emission lines with their ion 
species and redshifts in a similar fashion to Figure [S] 

In Figure [7|i it is clear that the IGM emission domi- 
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Fig. 6. — Composite spectrum of the Galactic foreground and 
AGNs. It includes three components: (a) a thermal emission from 
the Local Hot Bubble; (2) an absorbed thermal emission from hot 
gas either in the distance halo or in the Local Group; and (c) 
continuum emission from discrete AGNs. 



nates the entire spectrum. This is not surprising because 
the X-rays come from the hot gas within a galaxy group 
at redshift z w 0.07. For such hot gas, bremsstrahlung 
radiation, which dominates the total X-ray emissivity, is 
oc n^T^-^, where Ue is the electron density and T is the 
temperature. The X-ray emission from the high density 
and high temperature intragroup medium shows a strong 
continuum, and this dominates over the AGN and Galac- 
tic foreground emission. In this panel we identify emis- 
sion lines from O, Fe, and Ne. In particular we can easily 
make out redshifted O VIII Lya (resonance line at 611 
eV) and the O VII triplet (forbidden line at 524 eV, in- 
tercombination line at 531 eV, and resonance line at 536 
eV). This could provide us with a direct measurement of 
the temperature, electron density and ionization states 
of this X-ray emitting gas. The numerous lines between 
650 - 850 eV are mainly emission lines from Fe L-shell 
transitions, and many of these lines belong to Fe XVII 
transitions. 

Figure shows the emission spectrum of a patch of 
sky covered with more diffuse structures, known to be 
filamentary in 3 dimensions [area (b) in Figure The 
total spectrum (black line) is dominated by the back- 
ground (AGNs) and foreground emission, especially at E 
< 0.4 keV. However, we see several emission lines from a 
filament crossing the line of sight at z « 0.21. Although 
the overall continuum level is nearly an order of magni- 
tude lower than that of the AGNs plus the foreground, 
these emission lines are strong enough to poke through 
the total spectrum and show up as obvious features. 

In Figure [7|; the total spectrum is completely domi- 
nated by the AGNs plus the foreground. The emission 
from the IGM is almost 100 times lower that that of the 
AGNs plus the foreground, and it is not possible to detect 
any signals from the IGM. 

In Figure IHl and El we zoom in on areas (a) and (b), 
respectively, plotting only emission from gas which lies 
within particular narrow range in redshift (explained be- 
low). These plots have angular resolution four times 
higher than that in Figure El Figure El shows a sim- 
ulated 20' X 20' region, centered on area (a) (denoted 
by the white box). Most X-ray emission lines from this 



area originate from hot gas at z « 0.0775, so that we 
have cut out a slice that is centered on this redshift with 
Az = 0.0005 (about 1.5 comoving h'^Mpc). Figure El 
shows the X-ray emission from this slice, and the dotted 
lines show contours of X-ray emission. In similar fashion 
we have made Figure for area (b). Clearly, the round 
shape of the X-ray contours in FigureElindicates that the 
hot gas in area (a) sits in a virialized, relatively stable 
system, such as a galaxy group. The X-ray contours in 
area (b) (Figure EJ (we are plotting only emission com- 
ing from close to z = 0.21) show rather more diffuse and 
irregular X-ray emission from filamentary structures. 

To further understand the gas distribution within area 
(a) and (b), we examine all the particles in these two 
areas centered around redshifts of ~ 0.775 and 0.21, 
respectively. The length along the line of sight corre- 
spond to 5' at that redshift, i.e., two cubes with comov- 
ing sizes 0.33ft,-i x 0.33/1"! ^ 0.33"! Mpc^ at z 0.775 
and 0.86/i~i x OMh~^ x 0.86-iMpc3 at z = 0.21. In 
Figure irni we plot the temperature vs. overdensity of gas 
particles with T > 10^ K. Red filled circles represents 
particles in area (a), and black ones are from area (b). 
In area (a), most particles concentrate on a region with 
overdensities around 10'^, while in area (b) gas particles 
distribute more diffusely. The confirms our initial guess 
that particles in area (a) are more concentrated and tend 
to form a virialized structure, and in area (b) structure 
is more filamentary-like. 

Clearly, an important method to distinguish between 
area (a) and (b) is to study the emission line strength 
via, i.e., Figured But exactly how does the observed 
line strength vary with gas density? To understand this, 
we select O VIII line at 654 eV as a representative line 
and study the line flux in each of the 5' x 5' grid. Again we 
select all the grids around the two redshifts, 0.0775 and 
0.21, to calculate the O VIII flux from each cube. The 
total flux of O VIII emission line from one cube can be ob- 
tained by direct summarization of the line flux from each 
individual particle in that cube l|Yoshikawa et alJ l2003). 
The mean density of that cube is calculated by averag- 
ing the total mass of particles with T > 10^ K within 
that cube. Figure [TTl shows the line flux distribution for 
both regions: red flUed circles represent cubes at redshift 
of ^ 0.0775, and black filled circles are from cubes at 
z 0.21. Two circles pointed by arrows are area (a) 
and (b). While the line fiux dependency on gas density 
is a complicated relation because it involves properties 
like temperature, ionization structure and metallicity, in 
general Figure^Jshows that high density regions tend to 
have high O VIII fiuxes. We also calculate the flux sensi- 
tivity of the proposed X-ray mission MBE (the horizon- 
tal dashed line) for a 200 ksec exposure time. While only 
very few girds in this plot are detectable with MBE at 
the two redshifts, the cumulative effect by observing each 
grid up to higher redshifts greatly enhances the overall 
detectability (see discussion in section §5). 

4. CORRELATION 

4.1. Angular Clustering of the WHIM 

As most of the emission from the WHIM comes from 
relatively low redshifts, z < 0.3, the strong spatial clus- 
tering of the baryons on small scales is directly reflected 
in strong angular clustering of their X-ray emission. In 
Figure 1121 we show a map of a region of sky of angu- 
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Fig. 7. — Simulated spectra of three selected regions. Green lines are for the IGM emission, red lines are for the Galactic foreground plus 
AGNs, and black lines stand for total X-ray emission, (a) Hot gas in a galaxy group, it is clear that the IGM emission dominates the whole 
spectrum, especially at high energy; (b) diffuse X-ray emission in a filament; The Galactic foreground plus AGNs' emission dominate, but 
some strong emission lines from a filament at z fs; 0.21 show up in the total spectrum; and (c) a void, in which the entire spectrum is 
dominated by the Galactic foreground plus AGNs. We also label redshifts of several selected lines. 



lar size 10° x 10° taken from om' hydrodynamic sim- 
ulation. The image is taken in the energy band from 
0.5 - 0.8 keV, and includes an AGN contribution and 
a galactic foreground component (both smoothly dis- 
tributed). The mean flux from AGN plus foreground 
is ~ 4 X 10~^^ ergs cm~^s~^deg~^, which we obtained 
by integrating the model from Figure between 0.5 - 
0.8 keV. The map has been smoothed with a filter of 
Gaussian width a — 5'. In the map, emission from many 
groups of galaxies can be seen as hotspots, as well as 
more diffuse emission between them. 

The angular clustering of the WH I M em ission is 
treated more extensively in ICroft et alJ l|200H) . In this 
paper we show an angular map with a larger scale, and in 
a narrower band. It is evident that the WHIM is strongly 
clustered. We compute the angular correlation function 
w{Q) from this map, where w{&) is defined as: 



w{d) = (,5x(a)<5x(a + 0)) 



(2) 



Here Sxict) is the X-ray flux over the mean: &x = 
fx/ (fx) — 1 where fx is the X-ray flux in each pixel 
between 0.5 - 0.8 keV. We look at three types of X-ray 



emission auto-correlation function: X-rays from WHIM 
only, in which we solely use gas particles with tempera- 
tures between 10*" - 10^ K; X -rays from the emission lines 
only, for which we subtract the continuum from the total 
spectrum; and X-rays from the IGM. All three types ex- 
hibit a strong clustering signal at high significance (see 
Figure lT^ . This strong clustering is an obvious signature 
of the IGM. Photons from AGN originate from much far- 
ther away and their angular clu stering is predicte d to be 
smaller on these scales (see e.g. . ICroft et al.l[2001|) . Even 
with only one such field, the IGM signal shows up very 
well. 

4.2. Cross- correlation of the IGM Emission and Galaxy 

Positions 

Numerical simulations show that the X-ray emission 
from the IGM traces the same large scale structures as 
galaxies. By calculating the angular cross-correlation of 
the X-ray background (XRB) with nearby galaxies in 
simulated surveys of 1° x 1° fields, ICroft et all l|2001fl 
found a strong signal on an angular separations of < 5' 
(see their Figure 13a). Since the X-ray emissivity is 
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Fig. 8. — A slice of the simulated region, centered on area (a) 
(white box), with a size of 20' X 20'. The slice is a projection 
of Az = 0.0005, centered at z fa 0.0775. The X-ray contours are 
separated by a factor of ~ 2.2 in intensity, with the highest contour 
at a level of ~ 10"''^ ergs s~^cm~^deg~ . 
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Fig. 9. — Same as Figure ISl but enlarged area (b) (white box 
at the center). The slice is a projection of Az = 0.0005, cen- 
tered a.t z ^ 0.2108. The X-ray contours are separated by a fac- 
tor of ~ 3.5 in intensity, with the highest contour at a level of 
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oc ri^T"', and the temperature is roughly propo rtional 
to the density for WHIM gas ijPave et al.l l200H) . it is 
obvious that the large scale pattern of X-ray emission 
should correlate well with that of baryonic overdensity 
S. Moreover, based on the assumptions in our simula- 
tion, the highly ionized metals can be directly related to 
local matter overdensities, so that we expect a strong cor- 
relation between emission lines from metals and nearby 
large scale structures. In this section, we will test the 
idea that correlation of spectral pixel fluxes with the po- 
sitions of galaxies will enable us to detect the presence 
of the WHIM. 

We first identify galaxies in the simula- 
tion using the SKID groupfinder (see, e.g., 
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Fig. 10. — Temperature vs. overdensity of gas particles with 
T > 10^ K in area (a) and (b). Red filled circles represents particles 
in area (a), and black ones are from area (b). 
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Fig. 11. — O VIII line flux vs density in each cube at redshift 
0.0775 (red) and 0.21 (black) in our 200' X 200' simulation. We 
also label the two areas (a) and (b). The dashed horizontal line is 
the estimated flux sensitivity of the proposed X-ray mission MBE 
for a 200 ksec exposure time. 



iKatz. Weinberg. &: Hernouistl Il996j) ". The galax- 
ies are then used to create a mock redshift survey 
occupying the same simulation space as the X-ray 
emission. T his part of the pro cedure is similar to that 
employed in lCroft et al.l l)2001(l . and we also produce a 
flux limited survey of galaxies. We adjust the magnitude 
limit so that the mean redshift of the galaxies is z = 0.2. 
We apply an upper redshift cutoff to our sample at 
z < 0.5, because most X-ray emission from the hot IGM 
comes from low redshift. The simulated sky area which 
we explore in this section has a size of 200' x 200', and 
the redshift survey contains 57,607 galaxies. 

We divide the 200' x 200' field into a 200 x 200 grid, 
with each cell subtending 1' x 1'. We calculate the to- 
tal spectrum in each cell, accumulating emission from a 
redshift up to z ^ 0.5, and using a spectral resolution of 
2 eV. 

Our analysis technique involves first subtracting the 

see |http://www-hpcc.astro.washington.edu/tools/skid.html| 
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Fig. 12. — A simulation of 10° X 10° region. The image is an 
integration of the total spectrum (IGM plus Galactic foreground 
plus AGNs) between 0.5 - 0.8 keV. The map has been smoothed 
with a filter of Gaussian width cr = 5'. In the map, emission from 
many groups of galaxies can be seen as hotspots, as well as more 
diffuse emission between them. 
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Fig. 13 . — Angular auto-correlation function calculated from 
Figure 1121 X-ray emission has been integrated from 0.5 - 0.8 
keV. Open triangles represent X-ray emission from the WHIM gas 
(10^ < T < 10^ K), open circles are for emission lines only, and 
filled circles are for X-ray emission from all of the IGM. 



continuum from the spectrum to obtain an emission-line 
only spectrum. We do this using an iterative technique 
to lower the emission line profile to the level of the con- 
tinuum, without using knowledge of the true continuum. 
To study the cross-correlation between galaxies and the 
X-ray emission lines from different ion species, for each 
ion species at a time, we remap all the energy bins in 
the spectrum, which contains emission lines from all ion 



species, to redshifts, based on the rest energy of that par- 
ticular ion species. For instance, to examine the cross- 
correlation between the O VIII emission lines and nearby 
galaxies, we set the zero-redshift point of each spectrum, 
which contains emission lines from all ion species, to 
Eq = 654 eV. The redshift Zi of spectral bin with en- 
ergy Ei is assigned according to z = Eq/E.^ — 1. Then for 
each spectral bin with energy Ei, we have a correspond- 
ing redshift Zi and two-dimensional angular coordinate. 
We then use the redshifts of the spectral bins and their 
angular positions on the sky to give each spectral bin a 
three-dimensional comoving Cartesian coordinate (mak- 
ing use of the appropriate relations for our ACDM cos- 
mology.) We define the cross-correlation function £,c{R) 
to be: 

UR)^{Ix{r)Sg{r + R)). (3) 

Here Ix is the X-ray intensity in each spec- 
tral bin after continuum subtraction, in units 
of photons cm~^s~^sr~^keV^^, and 6g is the 
galaxy overdensity. S,ciR) therefore has units 
photons cm^^s^^sr^^keV^^. 

Figure El shows the cross-correlation functions for 
eight ion species. During the calculation, the zero- 
redshift points were set to the corresponding rest-frame 
energies, as mentioned above. For He-like ions, we chose 
the rest resonance energies, and for Fe XVII we let 
Eq = 825.79 eV. The 1-a error bars are obtained by 
splitting the dataset into 20 subsamples, and calculating 
the standard deviation of these 20 subsamples. It is ap- 
parent that for all the ions except C V and N VI, there 
are strong signals at small distances, which we interpret 
to mean that these X-ray emission lines are strongly cor- 
related with large scale structures traced by galaxies. In 
regions close to galaxies, these signals, after subtracting 
their mean values at large distances (by analogy with the 
galaxy or matter correlation function), can be fitted with 
power laws with slopes ranging between -1.2 - -2. Since 
hot galactic halo gas lies typically within scales less than 
a few hundred Kpc from galaxies, we can identify these 
signals as the correlation signals of hot gas in the IGM 
and galaxies. If we do not subtract the mean value of 
(^c) from the results, the log(^c) ~ log(i?) relation can 
be approximately fitted by a power law over the entire 
range. Fe XVII and O VIII show the steepest power law 
slopes of - (-0.43), while O VII, N VII, and C VI show 
slopes between (-0.33) - (-0.30). Ne IX has a rather flat 
slope of ~ (-0.15). 

It is interesting to see that different ions show different 
characteristic correlation strengths. The strongest signal 
comes from O VIII and Fe XVII. O VII and Ne IX also 
show a relatively strong signal, whereas C VI and N VII 
give a weak signal, while C V and N VI are not detected 
at all. This can be explained by the combination of ion 
emissivity and element abundance. From our definition, 
the correlation function £,c{R) should be a function of 
(eiZi), where ei is the line emissivity for ion Xi, and 
Zi is its elemental abundance. The higher (tiZi), the 
larger the £,c{R)- From Figure^ O VIII (solid blue-green 
curve) and Fe XVII (dotted green curve) show the highest 
peak emissivities, so they have the highest correlation 
signals. However, since oxygen is more abundant than 
iron, O VIII gives a stronger signal, even compared to 
Fe XVII. Following a similar line of argument, we can 
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Fig. 14. — Cross-correlation between the IGM emission lines and galaxies: (a) Ne IX and Fe XVII, (b) O VII and O VIII, 
N VII, (d) C V and C VI. These IGM emission lines (except N VI and C V) tend to cluster around galaxies at scales ^ 1 - 



(c) N VI and 
2 Mpc. 



explain the relative strengths of the correlation signal 
from the other ions. 

Different ions are detectable over different scales, de- 
pending on the strength of the emission. The correlation 
signal for O VIII can be seen strongly for i? < 2 Mpc, 
while others (namely O VII, Fe XVII, Ne IX and C VI) 
only start to show an obvious signal at i? < 1 Mpc, with 
N VII only tentatively seen, at an even smaller scales. 
For highly overdense regions, which typically also have 
large galaxy overdensities, even on relatively large scales, 
the high elemental abundance and large emissivity gives 
a significant signal for O VIII; while for other ions, the 
lower abundances and/or lower emissivities mean that 
the signal is only readily seen within short distances from 
the overdense area, say, i? < 1 Mpc. 

5. DETECTABILITY 

Generally, most emission lines from the filaments are 
weak, undetected by current X-ray telescopes. These 
lines are expected to be embedded in the strong contin- 
uum emission from AGNs plus the Galactic foreground 



(see Figure E}. To investigate their detectability, we de- 
fine the equivalent width of an emission line at energy E 
as 

EW= I l\dE, (4) 

Je~AE/2 / 

where / and Ic are the observed spectrum 
and its continuum, respectively, in units of 
photons cm~^s~^keV~^sr~'^. The integration is 
over a small energy range AE around E. For a weak 
emission line, EW <^ AE. Given an exposure time of 
T, the minimum detectable equivalent width is 

Here {S/N) is the desired signal-to-noise ratio, R is the 
instrumental resolving power, A is the effective area and 
Q is the solid angle over which the emission is integrated. 

It is apparent, from equation jSJl, that an instrument 
with high {RAfl) is needed to detect weak emission fea- 
tures. Current X-ray telescopes, such as the Chandra 
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X-ray Observatory and the X-ray Multi-Mirror Mis- 
sion (XMM-Newton) have unprecedented spectral 
and spatial resolution. However, the detectors on board 
these telescopes either have high spectral resolution but 
low effective area, or large effective area but low resolv- 
ing power, and thus are inadequate for achieving high 
resolution spectroscopy of extended structures that we 
are interested in. 

Several proposed and planned X-ray missions have a 
large combined {RAQ) and will show promise in de- 
tecting weak emission features from extended structures. 
We now examine four future X-ray missions, namely 
Astro-E2 Constellation-X l^ The X-Ray Evolving 
Universe Spectrometer (or XEUS ^^), and the Missing 
Baryon Explorer (or MBE i"^). These X -ray missions are 
either in the final stage (Astro-E2), in planning stage 
{Constellation-X,XEUS) , or in proposal stage (MBE). In 
particular, we will concentrate on MBE, a mission specif- 
ically designed to detect diffuse X-ray emission from the 
WHIM, to see how this can be achieved. 

Table 2 gives the relevant parameters for these four in- 
struments. In calculating the EW, we assume that the 
desired S/N ratio is 3, that the line has energy 0.6 keV 
(where most O and Fe lines lie), and that the exposure 
time is 200 ksec. Since the observed spectrum is dom- 
inated in most cases by the emission from the Galactic 
foreground and the AGN background, we adopt the con- 
tinuum value Ic ~ 24 photons cm~^ s~^sr~^keV~^ from 
Figure El We assume that we will devote a time, T, for 
each telescope to observe a region of sky of solid angle, il, 
and we calculate the minimum EW of detectable emis- 
sion lines integrated over a solid angle, ^pixel- When a 
telescope with a field of view, QpoVj maps a region with 
solid angle, Q, with a total exposure time, T, each point 
of the region receives a exposure time of (T x flpoy/fl), 
so the minimum EW in each pixel can be expressed as: 
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Fig. 15. — MBE view of the same area shown in Figure [T^ 
The total exposure time is ~4 Msec. A Poisson-distributed noise 
was added. AGNs plus the Galactic foreground contribute a mean 
count of ~ 80 counts pixel 



Tabic 2; Instrument Parameters and Detection Limits 



Instrument 


A 




f^PIXEL 


R 


RAQ rn-iFOV 


EW 




(cm^) 












(cm^deg^) 


(eV) 


Astro-E2 ^ 


35 


2.9' 


X 2.9' 


0.48 


X 0.48' 


100 


5 


138.0 


ConstcUation-X ^ 


3,000 


2.5' 


X 2.5' 


5' 


X 5" 


400 


2,083 


8.0 


XEUS ^ 


40,000 


1' 


X l' 


2' 


X 2" 


500 


5,556 


4.9 


MBE 


300 


29.5' 


X 29.5' 


4.9 


X 4.9' 


150 


10,878 


3.7 



1. XRS: tlie X-ray Spectrometer 

2. SXT; the Spectroscopy X-ray Telescope 

3. STJ I: the initial configuration for the Supercon- 
ducting Tunneling Junctions. 



EW > 



QE 



IcRATQfov^pixel 



(6) 



To compare the minimum detectable EW among the 
four instruments, we use as the value for pixel the 
largest instrument pixel size (4.9' x 4.9' from MBE), and 
for we take the largest instrument field of view (29.5' x 
29.5' from MBE). In the case of the same il, ^pixel, 
and T, equation (5) indicates that the instrument with 
the largest value of RA^l,pov gives the smallest EW, as 
indicated in the last two columns of Table 2. Although 
XEUS has much higher angular resolution than MBE, 
its smaller instrument pixel size offers no advantage in 
detecting faint diffuse emission lines, unless the relevant 
scale of angular structure is significantly smaller than the 
MBE pixel size. The large field of view of MBE makes it 
clearly the best-suited mission for detecting weak X-ray 
emission lines from extended sources. 

see'http:/ /asc. harvard. ed u| 

see http://xmm.vilspa.esa.es3 



see http: / / www. isas.ac.jp/ e/enterp / missions^ast^>eii 
see http://constellation.gsfc. nasa.gov/docs^namTrtnu 
see http://astro.estec.esa.nl/XEUS 
^"^ see http://www.ssec.wisc.edu/baryonsj 



MBE is a mission designed specifically to observe the 
X-ray emission from the warm- hot, moderate overdense 
gas predicted to be distributed in filamentary structures 
that connect collapsed, virialized regions such as groups 
and clusters of galaxies. It will also be well suited for 
detailed studies of emission lines from groups themselves. 
It consists of a single instrument, the X-ray Calorimeter 
Telescope (XCT) - a high-resolution spectrometer and a 
moderate resolution imager. The energy resolution is ~ 
4 eV (FWHM) , with an energy bandpass from 40 — 2000 
eV. The field of view consists of a 6 x 6 array of pixels 
~ 4.9' on a side, so that the field of view is 29.5' x 29.5'. 

In Figure 1151 we show a simulated MBE count map 
of the same 10° x 10° field shown in Figure [T21 The 
energy band plotted is between 500 - 800 eV. The to- 
tal exposure time required to build up such an image is 
'^4 Msec, consisting of 400 pointings of 10 ks each. We 
add photons from AGNs plus the Galactic foreground, 
with a mean count of ~ 80 counts pixel^^. We also 
add Poisson-distributed noise to the map. During such 
a short exposure time, the counts collected in each pixel 
range from ~ 80 - 90 counts pixel^^, for pixels which 
contain purely background/foreground emission, up to 
^ 350 counts pixel^^, arising from hot spots such as 
galaxy groups. 

To demonstrate how MBE can detect weak X-ray emis- 
sion features, we also present the instrumental-folded 
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Fig. 16. — Simulated spectra observed with MBE response functions. The three spectra correspond to the three regions labeled in 
Figure|2l respectively. The exposure time is 200 ksec for each spectrum. In panel (b), it is clear that MBE should be able to probe several 
weak X-ray emission lines from filaments. 



spectra of the three areas we discussed previously (see 
Figure El- To obtain these spectra, we used the soft- 
ware packa ge ISIS (Interactive Spe ctral Interpretation 
System, see lHouck fc Denicolall2000 D and apphed the 
instrumental response data and effective area data for 
MBE. The three panels in Figure El represent spectra 
for the three regions we have examined before. The red 
lines are the model spectra, and the black crosses rep- 
resent observations of these spectra. The errors bars on 
each point are based on l-cr Poisson statistics only. The 
binsize is 2 eV and the exposure time is 200 ksec. We 
label the ion species and their redshifts in red text. It 
is clear that many weak emission lines, particularly in 
Figure ITCb . where the emission is from a filament, can 
be detected without ambiguity. For instance, the emis- 
sion line from the z « 0.21 O VIII line contains nearly 
30 counts. 

An interesting question to ask is what fraction of the 
sky as probed by MBE can yield detectable weak emis- 
sion lines such as those in Figure llfib . given a certain 
amount of exposure time? This we have calculated by 
examining the 1,600 spectra from the 200' x 200' sky 

18 see [httgr77ggi:ce.mrt .edu/ ASC/ISIS/ 1 



simulation. For most WHIM emission, the continuum 
should be dominated by the AGN plus Galactic fore- 
ground emission. Taking oxygen lines as an example, as- 
suming most weak emission lines seen in Figure span 
only one bin, which has a binsize of 2 eV, to ensure a A-a 
detection in one pixel, the signal in the line, as a function 
of exposure time T, should be at least 

S" = 1.85 X 10"^ counts, (7) 

and the line peak intensity, also as a function of exposure 
time T, should be at least 

/pg^fe = 2 X 10'' T'i + 24 photons cm"^s~^keV"^sr"\ 

Assuming an exposure time for each pixel, we can 
then convert this line peak intensity to line flux sen- 
sitivity. Taking T = 200 ksec as an example, the 
flux sensitivity at the energy of O VIII would then be 
1.4 X 10^^°ergs cm^^s^^sr^^. This is shown as the hor- 
izontal dashed line in Figure ^2 Typically MBE would 
be able to probe regions with overdensities around 100. 
In the 200' x 200' sky simulation, we can then count the 
total number of angular grid cells for which the spectra 
show at least one line with peak intensity > Ipeak- The 
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Fig. 17. — The fraction of the sky probed by MBE with at least 
one detectable (4(t) weak emission Une, given the exposure time T. 
The errors are based 1-cr poisson fluctuation only. 



fraction of the sky is obtained by dividing this number 
by the total number of angular grid cells, which is 1,600 
in this case. Figure El plots that fraction of the sky as 
a function of exposure time. The errors on the points 
are based on l-cr poisson fluctuations only. We can see 
from this that for one 200 ksec observation with MBE, 
given the field-of-view of ~ 900 square arcmin, roughly 
^ 5% of the field-of-view, or 2 pixels will show detectable 
emission lines. 

6. DISCUSSION AND SUMMARY 

We have presented a detailed study of predictions for 
the soft X-ray emission from the Warm-Hot Intergalactic 
Medium. Our main conclusions can be summarized as 
follows: 

1. We have examined images of soft X-ray emis- 
sion from the hot IGM with a wide field of view 
(200' X 200'). Having spectral information for each 
pixel allowed us to view the sky in several differ- 
ent soft X-ray energy bands, where emission lines 
from different ion species dominate. The low en- 
ergy emission is markedly more diffuse than X-rays 
in the higher bands. 

2. We selected three representative simulated patches 
of sky to study with high resolution X-ray spectra: 
(a) a galaxy group; (b) a filament; (c) a void-like, 
underluminous region. In each case, we saw nu- 
merous emission lines from elements such as C, N, 
O, Fe, Ne, and identified their redshifts. 

3. By taking into account the background X-ray emis- 
sion from the AGNs and foreground emission from 
the Galaxy, we obtained composite X-ray spectra 
of the selected regions. We found that in areas 
where the galaxy group is present, the X-ray spec- 
trum is dominated by the hot intragroup medium, 
while in the void-like areas, the spectrum is com- 
pletely dominated by the AGNs plus the Galac- 
tic foreground. In the filament case, although the 
majority emission comes from the AGNs plus the 



Galactic foreground, it is still possible to detect 
strong emission lines from the WHIM. 

4. We found that the emission in the (0.5-0.8 keV) 
band (covering many of the emission lines) shows 
strong small scale angular clustering, as measured 
by the angular autocorrelation function. Also, 
spectral information for each pixel allowed us to 
measure three dimensional clustering. In partic- 
ular, we computed the cross-correlation between 
eight different ion species and the simulated galaxy 
distribution. We found that O VIII showed the 
strongest correlation signal, being measurable on 
scales as large as ~ 4 — ^h"^ Mpc, while low en- 
ergy species such as C V and N VI showed no cor- 
relation at all. This variation of the correlation 
strengths can be explained by (eiZi), the combina- 
tion of line emissivity and elemental abundance 
Zf. the larger (eiZi), the stronger the correlation 
signal. 

5. Finally, we studied the detectability of the WHIM 
with several proposed X-ray missions. By com- 
paring the effective area, resolving power and in- 
strumental field-of-view, we found that the Missing 
Baryon Explorer (MBE) shows the most promise 
for detecting weak emission lines from the WHIM 
gas. We also showed simulated MBE instrument- 
folded spectra of our three selected patches of sky. 

Accurate prediction of the soft X-ray emission depends 
on several important factors, including X-ray emissivity, 
ionization fraction, metal abundances, etc. The largest 
uncertainty in our predictions of the soft X-ray emission 
comes from the metal abundance in the IGM. Theoreti- 
cally, simulations are beginning to address the formation 
and distribution of metals in a self consistent way (see 
e.g., ICen fc OstntoHTflflflt ISnringel HTT^i ist 200l 
Aguirre et al. 2001a, b,c). The metallicity-density rela- 
tionship we adopted in this paper was motivated by this 
work, and future studies of UV and soft X-ray emission 
lines should examine the q uestion of m etal en richment in 
more detail (see, e.g.. iFurlanetto et al.ll2003(l . 

Observationally, metals are known to be distributed 
in a variety of systems with large dispersion in metal- 
licity Z. For example, the metallicity can be as 
high as IOZq in the central regions of s ome active 
galac ti c nuclei (AGNs) dM ushotzkv. D one fc Pounds I 
119931 iTripD. Lu fc Savage! .1997) . and as l ow as 
10~^ Z(r) in some halo stars in our galaxy l|Beerg| 
119991) . When we focus on the low redshift IGM, 
the main evidence for an enriched IGM comes from 
the observatio ns of the X-rav-ern i tting gas within 
galaxy clu sters (|Mitchell et al. II1975I: i Mushotzkv et aT 
^96; Mus hotzkv fc Lo ewenstein |1997^|Dupke fc White 
mOOi) . With ^.SC^. IMu shotzkv et al. I (1199(11) found 
that the metallicities in four clusters vary between 
0.3 - 0.5Zp,. Studying ~ 40 c lusters at z > 0.14, 
IMushotzkv fc Loewensteinl l)1997!) concluded that there 
is no metallicity evolution in clusters up to z 0.3. 
The observations of the Lya systems at low redshifts 
give m etallicities as high as ^ 0.05^0 (Barlow fc TvtleO 
[1^)98) indicating a metal enrichment in even lower density 
systems at low redshifts. The observations of damped 
Lya systems [N(H I) > 10^" cm~^] show a similar result 
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ijde la Varga et all l200Cfl . Due to the lack of observa- 
tions, the metalhcity of the "WHIM" is unknown. 

Another uncertainty comes from the sensitive depen- 
dence of the X-ray emissivity and ionization fraction 
on physical properties of the WHIM gas, such as its 
temperature and density. For instance, both the emis- 
sivity and ionization fraction peak over a rather nar- 
row range of temperatures, and drop rapidly at low 
and high temperatures (see Figure [J and F igure 2, 
3, & 4 of IChen. Weinberg. Katz. fc Pavel [20f)l . which 
means a small change in temperature will cause order-of- 
magnitude changes in both quantities. This will in turn 
cause a dramatic change in observables such as intensity 
and emission line ion species. Because of this the actual 
observed spectrum depends sensitively on the underlying 
physical properties of the WHIM gas. 

In this paper, we adopted several approximations to 
simply the calculations. First we did not take into ac- 
count photoionization - we used coUisional ionization 
when calculating the ionization fractions. It is impor- 
tant to be aware that photoionization from background 
radiation may substantially alter the metal ionization 
fractions for low density, low temper ature gas (see, e.g., 
iChen. Weinberg. Katz. fc Davell2003() , which in turn will 
change the observed spectrum. However, since most of 
the emission spectra that can be observed with MBE 
or other telescopes are produced by hot gas in the high 
density, high temperature tail of the WHIM distribu- 
tion, where coUisional ionization dominates, we expect 
that photoionization will not substantially change those 
observables. 

Metal cooling is ignored in our simulation: the cooling 
is primarily through H and He. Metal line cooling could 
be very important in some c ases, especially for high den - 
sity, warm/hot gas (see, e.g.. lSutherland fc DoDitall993j) . 
However, since our results are very sensitive to the tem- 
perature variation, it is important to quantify the possi- 
ble effect of metal cooling on WHIM gas temperature. To 
cool some hot gas, giving it a small drop in temperature 
ST, the characteristic time scale St is 



St = 



nA(T) 



ST. 



(9) 



Here a is 3/2 or 5/2, depending on whether the cooling 
is isochoric or isobaric; k is Boltzmann's constant; n is 
the electron density; and A(r) is the cooling function. 
Assuming that we are interested in how cooling acts on 
timescale of order the Hubble time, we can estimate the 
effect of metal cooling on the temperature decrement ST 
in terms of the physical parameters of the WHIM gas: 

ubblc; 

(10) 

where n_5 = n/(10-^ cm-^); Tg = r/(106 K); Stanhhie 
is the cooling time in units of Hubble time scale. Here 
A_23 is cooling rate in units of 10~^^ ergs cm^^s"^, 
typical of the WHIM gas with 0.1 solar abundance 
ijSutherland fc Dopita"1993'). So in general, metal cool- 
ing has only a minor effect on our results; however, since 
A(T) oc T^^ in the temperature range we are interested 
in, the metal cooling process is more efficient for low tem- 
perature gas, especially so when A(r) reaches its peak 
at T w 2 X 10^ K. 

It has been known for quite some time that non- 
gravitational heating, and/or cooling processes can 



change the ov erall l evel of cosmological X-ray emission 
substantially llPenI Il999t \Wu. Fabian, fc NiilsenI 120011 
IVoit fc Brvanll2001aD . The energy deposited by super- 
novae in our simulations has little direct effect on the 
IGM. It is added as thermal energy to particles in high 
density regions and is radiated away quickly. Although 
the direct impact of feedback is therefore minimal, the 
cooling included serves to change the entropy distribu- 
tion of gas near the centers of virialized objects and 
lower the overall level of emission. Other treatments 
of starburst energy such as those which involve giving 
wind particles kinetic energy so that they escape be- 
yond galaxy and gro up vir ial radii into the IGM (e.g., 
Springcl fc Hcrnauis3 l2003j) are also likely to change the 
distribution of X-ray emitting gas. Different treatments 
of feedback have been shown by Brvan fc Vo it ( 2001) and 
iVoit fc BrvanI l)2001bj) to change the pdf of X-ray surface 
brightnesses in the range which would affect the WHIM 
emission lines we focus on here. Varying the feedback 
model would be useful in the future in order to assess 
the range of effects on the WHIM emission. 

Finally, the ionization state of the gas was calculated 
assuming ionization equilibrium. Non-equilibrium evolu- 
tion of the ionization state ca n significantly change the 
ionization fraction (see, e.g., lEdgar fc ChevaliCT Il986t 
Hicks fc Canizarcs 2OOI4). The characteristic time scale 
for hot gas to reach ionization equilibrium, or the re- 
combination time scale is ~ {narec)~^, where arec is 
the recombination rate. At a temperature of ~ 10^ 
K, the recombination ra te for O V II, for example, is 
- 1.4 X 10^12 cm^s-i l|Shull fc van Stccnbcrg 1982). 
Here we include contributions from both radiative and 
dielectronic recombination. Assuming n ~ 10~^ cm""^, 
typical of WHIM gas, the recombination time scale is ~ 
2 Gyrs, less than the Hubble time. However, we must 
keep in mind that most of the hot gas was shock-heated 
at low redshift, so it is very likely that a significant por- 
tion of the WHIM gas has not reached its recombination 
time scale and is in a non-equilibrium ionization state. 

A reasonable question to ask is: given detected X- 
ray emission lines, how do we identify lines with ion 
species and redshifts? It is impossible to carry out 
line identification with a single line; however, by simul- 
taneously observing multiple lines such as the O VII 
triplet and/or O VIII redshifted to the same degree, 
which is within the capability of the energy resolution 
of proposed X-ray missions such as MBE, we can eas- 
ily identify lines with redshifts. In addition, measure- 
ments of the line ratios can provide direct tempera- 
ture/density diagno stics for t he ionized plasma (see, e.g., 
Veddcr. Canizarcs. Markert. fc Pradhan 1986). Obser- 
vations of soft X-ray emission lines have the potential 
to not only reveal the presence of the missing baryons 
but to give us detailed information on their clustering, 
temperature, and metallicity. 

X-ray emission provides a direct probe of the WHIM 
gas. Unlike the X-ray absorption method, which can 
yield only one dimensional information, by combining 
imaging/spectroscopic observations. X-ray emission can 
show us the fully three-dimension distribution of the 
WHIM gas. This is essential not only for closure of the 
local cosmic baryon budget, but also if we are to un- 
derstand the distribution of the WHIM gas, and how 
structures form and evolve. Moreover, by applying line 
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diagnostic methods, X-ray emission lines can be used to 
probe the metaUicity, temperature and density of X-ray 
emitting gas. This in turn can give us important infor- 
mation on galaxy formation and evolution, given that 
most of these metals which are responsible for X-ray 
emission lines are produced by supernova explosions and 
propagate from the interstellar medium into intergalac- 
tic medium. Although the relatively low density and 
temperature make it hard to observe, currently proposed 
X-ray missions will for the first time provide us with de- 



tectable signals from the X-ray emission of the WHIM 
gas. 
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